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Abstract 

Recent results from the HI and ZEUS Collabora- 
tions on inclusive single and multiple jet production 
in the neutral current deep inelastic scattering of 
electrons/positrons on protons and the high energy 
photoproduction on protons at HERA are reported. 
The results are compared with NLO QCD calcula- 
tions and have been used to determine the strong 
interaction coupling constant as- 

1 Introduction 

The study of jets in the deep inelastic scattering 
(DIS) of electrons/positrons on protons and in high 
energy photoproduction on protons enables sensi- 
tive tests of QCD to NLO and provides precision 
measurements of the strong interaction coupling 
constant ag. In this report a survey is given of 
results from recent measurements of jet production 
in neutral current processes from the HI and the 
ZEUS Collaborations at the HERA collider; not in- 
cluded in this review are results on jets associated 
with diffraction, heavy flavours or direct photons. 




Figure 1: Examples of leading-order contributions to 
inclusive jet production: (left) direct photon interac- 
tion, (right) resolved photon interaction. 



In jet studies QCD predictions for DIS and pho- 
toproduction are probed in different regions of the 
DIS variables (Bjorken)x and (with y — / sx 
and s the squared ep centre-of-mass energy) and 
the transverse energy and pseudorapidity r]j 

*Talk given on behalf of the HI- and ZEUS- 
CoUaborations at the conference "New Trends in High- 
Energy Physics", Alushta, Crimea, Ukraine, May 24-31, 
2003 



(defined as rjj = — logtg 6j/2 with Qj the polar 
jet angle with respect to the proton beam) of the 
hadronic jet or jets. In DIS at high the virtual 
photon emitted from the electron interacts with a 
parton from the proton in a direct (pointlike) pro- 
cess while at low and in photoproduction it can 
fluctuate into partons, one of which interacts in so- 
called resolved processes with a parton of the pro- 
ton (fig.P). 

In pQCD jet production depends on the hard 
partonic scattering process and on the parton den- 
sity functions (PDFs) of the proton and, in photo- 
production and low processes, also of the pho- 
ton. For the calculation of the jet cross section aj 
in DIS (for photoproduction see Chap. 3), the hard 
partonic cross section daa for a parton a has to be 
convoluted with the PDFs fa of the proton 

= ^ / dxfa{x,^i%)daa{x,as{^J.R),^J.%^i%) 

a •' 

•(l+4adr) (1) 

and corrected for hadronization effects. The de- 
pendence of the PDFs on the characteristic en- 
ergy scale of the interaction is described by evo- 
lution equations i.e. for not too small x and 
by the DGLAP equations. The summed higher or- 
der contributions of this evolution, corresponding 
to multiple coUinear gluon emissions of the inter- 
acting parton, are represented by parton cascade 
diagrams. The factorization theorem allows the 
soft (nonperturbative) part of these contributions 
to be absorbed in the PDFs while the hard part 
has to be accounted for in the perturbative part of 
the interaction, it defines the factorization scale fip- 
Examples of LO contributions to DIS are shown in 
fig. 12 A breakdown of the DGLAP approxima- 
tion is to be expected in the small x region, where 
log l/x-terms become larger than log Q^-terms; in 
this region, the BFKL evolution equation should 
give a better description of DIS. A kind of inter- 
polation between these approximations is provided 
by the CCFM approach. 

A preferred reference system for analyzing jet 
processes in DIS is the Breit system in which the 
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parton collides head-on with the timelike virtual 
photon: 2xp+q= 0. In the quark parton model 
(QPM) the hadron jet resulting from the backscat- 
tered quark emerges collinear with the proton rem- 
nant, whereas QCD hard processes are character- 
ized by jets with high transverse momenta with 
respect to the proton remnants. Analyzing the 
hadronic final state in the Breit system therefore 
suppresses the QPM background. 

The identification of jets is achieved by algo- 
rithms which are applied to the objects of the mea- 
surement i.e. the particle tracks and energy clusters 
that are registered by the detector. Different jet al- 
gorithms have been developed and applied to the 
HERA experiments but in recent years the "inclu- 
sive longitudinally invariant fcx cluster algorithm" 
(in the following: ILICA) has become the preferred 
one. The algorithm proceeds by calculating for 
each track or energy cluster i the quantity di and 
for each pair the quantity dij : 

di — -E'T,i > 

d,, = mm{Eh,El^){A7^l+Acl,l) (2) 

where -Et,! is the transverse energy of a particle i 
and Arjij and are the differences in the pseu- 
dorapidity r/ and azimuthal angle of i,j respec- 
tively; all these quantities are invariant under lon- 
gitudinal boosts i.e. boosts in the direction of the 
proton beam (defined as -l-z-direction). If of all re- 
sulting {di,dij} the smallest one belongs to {di} 
it is kept as a jet and not treated further if it be- 
longs to {dij} the corresponding tracks or clusters 
are combined to a single object. The procedure is 
repeated until all tracks or clusters are assigned to 
jets. The algorithm is infrared and collinear safe, 
and can be applied to measurements of DIS and 
photoproduction processes. It typically requires 
smaller hadronization corrections than other clus- 
tering algorithms. 

For tests of QCD, the measured jet cross sections 
have to be compared with calculations performed 
(at least) to NLO. These partonic final states have 
to be combined with models for the development of 
the parton cascade and the hadronization. Models 
for the parton cascade are the parton shower model, 
implemented in HERWIG and LEPTO, and the 
dipole cascade model, implemented in ARIADNE; 
the hadronization of the partons is modelled by the 
cluster model, implemented in HERWIG, or by the 
string fragmentation model, implemented in JET- 
SET and used in LEPTO and ARIADNE. In RAP- 
GAP the LO QCD matrix elements are matched 
to DGLAP based parton showers in the LL ap- 
proximation; it allows RAPGAP, besides simulat- 
ing direct processes, also the simulation of resolved 
photon processes. In CASCADE, the QCD matrix 
elements are combined with parton emissions de- 
scribed by the CCFM evolution equation with the 




Figure 2: Leading order diagrams for (di)jet produc- 
tion: (a) photon-gluon fusion; (b) QCD-Compton pro- 
cess; (c) parton cascade. 



unintegrated gluon density distribution used as in- 
put. 

In photoproduction, direct and resolved inter- 
actions are simulated with PYTHIA, HERWIG 
and PHOJET. The hard partonic interaction is 
described by LO QCD matrix elements and the 
parton cascade is simulated by initial and final 
state parton showers in the LL approximation. In 
PYTHIA and PHOJET, the hadronization is sim- 
ulated with the Lund string model as implemented 
in JETSET; in HERWIG the cluster model is used. 

Commonly used parametrizations for the PDFs 
of the proton are those of CTEQ and MRST, while 
for the photon PDFs the GRV and AFG sets are 
used. 

2 Jets in deep inelastic scat- 
tering 

Results on inclusive jet production in DIS have 
been obtained on single jets, dijets and trijets and 
recently also on subjets. Since jet production 
at small x in the forward direction is related to 
forward single hadron production, results on 7r°- 
production are also discussed. 

2.1 Inclusive jets 

At high i.e. for 150 GeV^ < < 5000 GeV^, 
the inclusive jet cross section has been measured 
by HI in the Breit frame as a function of the 
jet transverse energy Et for 7 GeV < Et < 50 
GeV, 0.2 < y < 0.6 and -1 < -quh < 2.5. Jets were 
identified using the ILICA. Over the whole range 
of Et and Q^, the NLO calculation corrected for 
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hadronization effects (< 10%) gives a good descrip- 
tion of the data (fig. Ell- 
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Figure 3: Inclusive jet cross section in DIS as fct. of 
tlie transverse jet energy in the Breit system Etje- 



Fitting the QCD prediction to these cross sec- 
tions using CTEQ5M1, fifj = Et and — 
V 200 GeV^ (the average Et of the jet sample) 
yields as- First in separate fits to the data points 
in the four regions ^^(E't) was determined and 
evolved to as{Mz)- The consistency of these fits 
justified a combined fit to all data points with 
the final result, taking into account the uncertain- 
ties of the renormalization scale {fiR — Exj) and 
from the parton parametrization, of as{Mz) = 
0.1186 ± 0.0030(exp) + 0.0039/ - 0.0045(theor) + 
0.0033/ - 0.0023 (PDF). The experimental error is 
dominated by the hadronic energy scale of the LAr 
calorimeter with the statistical error contributing 
only 0.0007; the theoretical uncertainty includes 
about equal contributions from the hadronization 
and renormalization scale uncertainty. The result 
has been shown to be stable against variations of 
the jet algorithm. 

Similar measurements have been performed re- 
cently by ZEUS For > 125 GeY^,ETB > 8 
GeV, —2 < TjjB < 1.8 and demanding for the angle 
7 of the hadronic system — 0.7 < cos..^ < 0.5, the 
Q^- and i?TB-dependence of the inclusive jet cross 
section is found to be in good agreement with NLO 
QCD calculations. 

For the determination of as, the differential 
cross sections da/dV, where V = Q'^^Etbj, were 
calculated in NLO QCD for three MRST99 sets 



i.e. central, as tt as H, taking for the 

partonic cross sections the as{Mz) of the corre- 
sponding PDFs. For each bin i in V, the calcu- 
lated cross sections were parametrized according to 
[daiasiMz))/dV]^ = Ci.y as{Mz) + C2,r al{Mz). 
The best combined fit {y^ = 2.1 for 4 data points) 
was achieved for > 500 GeV^ resulting in 
as{Mz) = 0.1212 ± 0.0017 (stat) 
data. 
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Figure 4: Inclusive jet cross section in DIS as fct. of 
the transverse jet energy in the Breit system Etje- 



The energy scale dependence of as has 
been investigated using the same procedure but 
parametrizing the measured cross section in terms 
of as{{E!^, ^)), where (E'r^. ^) ist the mean value of 
E^.^^ in each bin; = i?x-^^ was chosen here since 
it provides the better description of these data. 
Good agreement with the predicted running of as 
is found over a large range of E'r^. ^ . 

These data have also been used to study the az- 
imuthal distribution of jets in the Breit frame 
At LO, the dependence of the cross section on the 
angle 0^ between the lepton scattering plane and 
jets produced with high E^ takes the form 



da 



A + Bcos4>f +C cos 24>f 



(3) 



where the coefficients A,B,C result from the con- 
volution of the matrix elements for the partonic 
processes with the PDFs of the proton. The 
cos 0-term arises from the interference between the 
transverse and longitudinal components of the ex- 
changed photon (for <C J) whereas the cos 2ip- 
term results from the interference of the -I-1/-1 
-helicity amplitudes of its transversely polarized 
part; nonperturbative contributions arising from 
intrinsic transverse momenta of the partons have 
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shown to be small and are negligible at high E^. 
The measured normalized azimuthal distribution 
shown in fig.|Slis well reproduced by NLO calcula- 
tions with either fiR — E'j^ or Q. 
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Figure 5: Normalized differential cross section for in- 
clusive jet production in DIS with Etbj > 8 GeV and 
— 2 < TjB < 1.8; inner errors for statistical, outer errors 
for combined statistical and systematic errors. 



HI has extended the inclusive jet measurements 
to low g2 values of 5 GeV^ < < 100 GeV^ for 
0.2 < y < 0.6, Et > 5 GeV and -1 < r/i^b < 2.8 
The differential cross section as a function of 
Et is shown in fig. IHl for three bins of rj. The 
good agreement with NLO QCD predictions (for 
fiR = Et) found for the backward and central 
region is getting worse in the forward region i.e. 
of 7/ > 1.5, and for E-r < 20 GeV, where NLO 
corrections and uncertainties are becoming larger. 
For most of the data the theoretical uncertainties 
(given by the renormalization scale uncertainty) are 
larger that the experimental errors. Investigating 
the dependence of this discrepancy, it was found 
to originate from the < 20 GeV^ region, indi- 
cating that NLO QCD works reasonably well even 
in the forward region Tyiab > 1.5 as long as both E^ 
and are not too small. 

In the forward direction, jet production is ex- 
pected to be sensitive to the dynamics of the par- 
ton cascade at small x. As suggested by Mueller- 
Navelet, the contribution of DGLAP evolution to 
DIS should be suppressed by demanding of the for- 
ward jet E^ « while that of BFKL evolution 
should be enhanced by keeping large as fea- 

sible and x/xj small. 

HI has taken data on inclusive jet production in 
the forward direction [5] i.e. 7° < Qjiah < 20° with 
0.5 < Q^/E^ < 2 and Xj > 0.035 {Et > 3.5 GeV). 
The measurement and jet search with the ILICA is 
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Figure 6: Inclusive jet cross section da/dET in DIS in- 
tegrated over 5 GeV^ < < 100 GeV^ and 0.2 < 
Viab < 0.6 as function of the transverse jet energy 
Et in the Breit system for different ranges of rji^h'- 
forward (proton direction): 1.5 < rj < 2.8; central: 
0.5 < < 1.5; backward: -1.0 < 77 < 0.5. NLO 
QCD model: DISENT with CTEQ5M (solid line) and 
no hadronization corrections; hatched band: change of 
the renormalization scale by 1/2 and 2. 



performed in the lab frame with cuts 5 < < 75 
GeV^ and 0.1 < y < 0.7. The data (fig. [U are 
compared with four different Monte Carlos. The 
differential cross sections in x and Xj are up to a 
factor of two larger than the DGLAP prediction 
for direct photons only (RG (DIR)), but they are 
reasonably well described if a resolved photon is 
included. The description provided by the CDM 
model is also good while CCFM predicts too large 
cross sections. 

In forward 7r° production, the parton is tagged 
by a single energetic fragmentation product. While 
smaller forward angles than in jet production can 
be reached, the cross sections are lower and the 
hadronization uncertainties higher. 

Inclusive 7r°-production in DIS has been mea- 
sured by HI [H], with the following cuts: 0.1 < 
y < 0.6, 4.10-5 < a; < 6.10-3, 2 GeV^ < < 70 
GeV^ and tt° cms-momenta > 2.5 GeV, polar an- 
gles of 5° < < 25° (corresponding to the central 
region in the hadronic cms) and > 0.01; no ex- 
plicit cut on p^^o/Q^ was apphed. In fig. |H1 the 
differential cross section is compared with QCD- 
based models. The DGLAP prediction with direct 
photon interaction only is too low, but the inclusion 
of resolved photon interaction - albeit with a choice 
of the renormalization and factorization scale dif- 
ferent {Q^ + 4p^) than in forward jet production 
- yields a good description. BFKL based calcula- 
tion also describes the data reasonably well, while 
CCFM evolution fails at small x. 
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Figure 7: Inclusive forward jet cross sections in DIS 
for ptj > 3.5 GeV; 7° < Qj < 20°; Xj > 0.035 at 
tlie iiadron level. Models: CDM (Color-Dipole-Model: 
DJANGO + ARIADNE), RG (RAPGAP: DIRECT 
and DIRECT + RESOLVED), CASCADE (CCFM). 
Data are preliminary. 
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Figure 8: Inclusive forward n° cross section in 
DIS for transverse momenta in the 7-proton cms 
Pt^o > 3.5 GeV and 5° < Qtv° < 25° for different 
ranges of Models: DIR (RAPGAP for DIRECT 
proc), DIR + RES (RAPGAP for DIRECT + RE- 
SOLVED proc), CCFM (CASCADE), CDM (Color- 
Dipole-Model: ARIADNE), mod. LO BFKL (modified 
LO BFKL evolution). 



In a recent measurement by ZEUS [7] of inclusive 
jet production in DIS for > 25 GeV^, y > 0.4 
and with Etj > 6 GeV, —1 < rjj < 3, deviations 
of the data from NLO QCD predictions calculated 
in the lab frame to 0{al) where found to increase 
towards the forward direction (fig. IHl upper part) 
and small E^j, and x. In a further analysis of 
these data, the Muellcr-Navelet suggestion was ap- 
plied by combining the detection of a hard forward 
jet {Etj > 6 GeV, < 77^ < 3) with requiring the 
hadronic angle 7hadn corresponding to the direction 
of the scattered quark in the QPM, to be recon- 
structed in the backward direction of the detector 
(cos7hadr < 0) and thus suppressing QPM contri- 
butions to the forward jets. Comparing the cross 
sections for such defined jet configurations ("dijet 
phase space") with NLO QCD predictions calcu- 
lated to 0{al), better agreement (fig. El lower part) 
is found in the forward direction at the expense of 
a considerably larger uncertainty in the rcnormal- 
ization scale which swamps a possible BFKL signal 
in this part of the phase space. 

2.2 Dijets 

The LO QCD processes which contribute to di- 
jet production in DIS are QCD Compton scatter- 
ing (QCDC) and boson-gluon- fusion (BGF). In the 
high region, where the QCDC process is domi- 
nant and the PDFs are well constrained by fits to 
inclusive DIS data, dijet measurements with as as 
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Figure 9: Upper figure: Inclusive jet cross sections in 
DIS for Q2 > 25 GeV^ y > 0.04, Etj > 6 GeV and 
— 1 < 77j < 3 in comparison to 0(as )-NLO QCD cal- 
culations in the lab frame ("inclusive phase space", see 
text). Lower figure: Jet cross sections for the same 
kinematical conditions but constraining the jet topol- 
ogy to r]j > and cos^hadr < in comparison with 
0(as)-NLO QCD calculations in the lab frame ("dijet 
phase space"). 



input enable tests of pQCD; at low Q^, where BGF 
dominates, the comparison of measured dijet cross 
sections with QCD predictions for different PDFs 
has been used to investigate the gluon density dis- 
tribution. 

In dijet measurements by HI jlj and ZEUS |S] 
asymmetric cuts in have been applied to the 
two jets highest Et in the Breit system in order 
to avoid regions of the dijet phase space which 
are sensitive to soft gluon radiation. The data 
(fig. Uni are compared with DISENT NLO calcu- 
lations, which have been checked with MEPJET. 
For 0(10) GeV2 < < 0(10.000) GeV^ and 
Etb > 5 GeV, i.e. where NLO and hadronization 
corrections are small, the studied jet observables 
are reproduced within about 10%; for < 10 
GeV^, however, where NLO corrections are be- 
coming large, significant disagreement is observed. 
(For new data on multijct production by ZEUS see 
Chap. 2.3.) 

From the dijet fraction i?2-i-i(Q^) measured by 
ZEUS, as was determined using the same method 
as applied to inclusive jets (Chap. 2.1) with 
the result a^iMz) = 0.1166 ± 0.0019(stat) + 
0.0024/0.0033(exp) + 0.0057/ -0.0044(theor). Also 
the dependence of as on the scale Q is in good 
agreement with expectations. 
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Figure 10: (left) Inclusive jet (crtot) and dijet (cr2+i) 
cross sections in DIS, hadroniz. corr. C refers to dijets; 
(right): dijet fraction -R2+1, shaded band gives uncer- 
tainty of abs. energy scale of jets. 



Recently, low x dijet production has been stud- 
ied in more detail by HI '5' in the kinematic region 
5 GeV2 < Q2 < 100 GeV2, 10"^ < x < 10^^ 
and 0.1 < y < 0.7. For jets reconstructed in the 
hadronic cms with cuts of E^i > 7 GeV, Et2 > 5 
GeV within —1 < 77 < 2.5, the triple differential 
dijet cross sections in {x , Q"^ , ETmax) (fig- EJ and 
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(a;, Q^, Arj) do not show significant deviations from 
NLO calculations with nn = {Eti +£'t2)/2, hf = 
= vTOGeV^, in contrast to the ratio of the 
number of dijet events with azimuthal angle sepa- 
ration of less than 120° to all dijet events (fig. El- 
Better agreement for part of the phase space can be 
achieved by combining LO ME with direct photon 
processes and kx ordered parton emission. 




Figure 11: Inclusive dijet cross section in DIS for 
as a function of the maximum 
transverse jet energy -E-rmax in the hadronic cms com- 
pared to NLO QCD predictions (for CTEQ5M) cor- 
rected for hadronization. Outer light error band in- 
cludes hadronization (dark band) and renorm. scale 
uncertainties. 



2.3 Trijets 

While the cross sections for inclusive jet and dijet 
production in LO pQCD depend on 0(q!s) only, 
the trijet cross section already in LO is sensitive to 
0{al). 

Trijet production in DIS has been measured by 
HI [El for 5 GeV2 < < 5000 GeV^, selecting 
jets with Et > 5 GeV in the Breit frame. The mea- 
sured cross sections (fig. I13|l were compared with 
LO and NLO calculations using NLO JET with /i^ 
and fip set to the average transverse energy Et of 
the 3 jets in the Breit frame; for the proton PDFs 
CTEQ5M1 was taken with as{Mz) = 0.118. The 
inclusive trijet cross sections for invariant masses 
> 25 GeV as well as the trijet/dijet ratio, in which 
many theoretical and experimental uncertainties 
cancel out, agree well with the NLO predictions in 
spite of the large NLO corrections. Also the angu- 



lar distribution of trijets agrees well with the pQCD 
expectation. 
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Figure 12: Ratio S of number of dijet events with 
small azimuthal jet separation ((ji < 2n/3) to the to- 
tal number of inclusive dijet events as a function of 
(Bjorken)a:: compared to NLO QCD predictions and 
with error bands as in fig. 1121 



Recent studies of dijet and trijet production by 
ZEUS JX for 10 GeV^ < < 5000 GeV^ in the 
kinematic range 0.04 < y < 0.6, cos7hadr < 0.7, 
with cuts in the jet phase space of — 1 < 77iab < 2.5, 
Etb > 5 GeV for invariant masses of M2j^3j > 25 
GeV are shown in figs. 1141 andll5l Good agreement 
is found with the pQCD predictions after correcting 
for hadronization, i.e. to 0{al). 

2.4 Subjet multiplicities 

Reapplication of the jet algorithm with smaller res- 
olution scale Ucut to identified jets allows the res- 
olution of clusters within the jets which are called 
subjets. The development of the subjet multiplic- 
ity with Ucut can be calculated in pQCD for high 
i?Tj and not too low j/cut values, i.e. where frag- 
mentation effects are estimated to be small in NLO. 
Since the subjet multiplicity is mainly determined 
by QCD radiation in the final state, it only weakly 
depends on the PDFs of the proton. 

The mean subjet multiplicity has been measured 
by ZEUS [El in NC DIS processes with > 125 
GeV^ for Et > 15 GeV and -1 < 77^ < 2. The 
ILICA was applied to jets in the lab frame, where 
calculations of up to 0(q;^) can be performed. The 
mean subjet multiplicity < risbj > was determined 
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Figure 13: Inclusive trijet production in DIS for 5 
GeV^ < 0^ < 5000 GeV^; -1 < r]ub < 2.5; Mg^ > 24 
GeV. QCD: NLOJET. 
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Figure 14: Inclusive dijet and trijet cross section in 
DIS. The light shaded band indicates the calorimeter 
energy scale uncertainty, the hatched band the renor- 
malization scale uncertainty. 



for 5.10 ^ < Hcut < 0.1. The data are in good 
agreement with NLO calculations with fii^ = fip = 
Q using the CTEQ4M PDFs. (fig.|Il|). 

Using the method described earlier (see Chap. 
2.1), as was determined from these data for 25 
GeV< Et < 71 GeV and ycut = 10"^ yield- 
ing as{Mz) = 0.1187 ± 0.0017(stat) + 0.0024/ - 
0.0009(syst) + 0.0093/ - 0.0076(theor). While the 
experimental uncertainty is comparable to that re- 
sulting from jet measurements, the theoretical un- 
certainty is larger and dominated by NNLO terms. 

3 Jets in high energy photo- 
production 

In LO QCD direct as well as resolved processes (c.f. 
Chap. 1) have to be considered in the calculation 
of the jet cross section. If Xj is the fractional mo- 
mentum of parton a of the photon and Xp the same 
of parton b of the proton, the jet cross section is 
obtained from the convolution of the photon PDF 
/-yb(a;-y, /i^) and the proton PDF fpaixp, fi'^) with 
the hard partonic cross section d(Ja,b{x-y, Xp, a^, 

(Tj — ^ ^ J" J" dX'ydXpfpa{Xp^ fj, ) 
a.b 

■dcra,b{Xp, X^,as, ^^^) ■ (1 + (5hadr) (4) 

summed over all partons a, b from the photon and 
proton. For direct photon interactions /-y is a 5- 
function at a;.^ = 1. Since the transverse energy i?T 



of the jet or jet system is the only scale parameter 
available in the photoproduction of jets from light 
quarks, /i^ and are set to /i/? — ^ ^ — E^. 
In resolved processes photon fragments may inter- 
act with proton remnants causing a background 
of multiparton interactions ("underlying events"), 
which has to be corrected for; it remains below 10% 
for Et > 20 GeV. 

At HERA, the proton PDFs are, with the cur- 
rent statistics, predominantly probed in the range 
0.05 < Xp < 0.5, where they are well constrained 
by measurements of Fl". For the photon PDFs the 
probed range is 0.1 < < 1 in which quarks are 
well constrained for x-y < 0.5 by measurements of 
F2 at LEP, while the gluon is only poorly con- 
strained. Since jet production at HERA is sensi- 
tive to the gluon component in the photon at LO, 
it should allow to extend the study of the photon 
structure to higher scales. 

In single jet production, in comparison to dijet 
production, an increased kinematic range is accessi- 
ble, higher statistical accuracy can be achieved and 
infrared dangerous regions of phase space are not 
existent; on the other hand single jet photoproduc- 
tion is less sensitive to details of the hard scattering 
process. 

3.1 Inclusive jets 

HI and ZEUS have performed similar mea- 
surements and analyses of inclusive jet cross sec- 
tions for photons of Q'^ < 1 GeV^ within — iKrjj < 
2.5. The HI data for Et > 21 GeV cover a range of 
the hadronic cms-energy W^p of 95 GeV < W^p < 
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Figure 15: Inclusive trijet cross section in DIS; error 
bands as in fig. 1141 
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285 GeV. The measurements were extended down 
to Et > 5 GeV with a dedicated trigger for pho- 
tons of < 0.01 GeV^ in the restricted kinematic 
range of 164 GeV < W^p < 242 GeV. The ZEUS 
data were taken for 142 GcV < W-yp < 293 GeV for 
Et: > 17 GeV (see fig. I17|l . HI compares the data 
with the NLO QCD calculations of Frixione and 
Ridolfi, using CTEQ5M1 for the proton PDFs and 
GRV for the photon PDFs and setting fi^ — fip to 
the average transverse energy of the two outgoing 
partons; ZEUS compares with Klasen, Kleinwort, 
Kramer, using MRST99 for the proton and GRV for 
the photon with fin = fip = Etj. In both analy- 
ses good agreement with NLO QCD over six orders 
of magnitude of the cross section is found; the HI 
data for < q.OI GeV and 5 GeV < < 12 
GeV may indicate a cross section rising above pre- 
diction with increasing rj. 

From the ZEUS data, as{Mz) was determined 
(see Chap. 2.1) by calculating in NLO QCD for 
each E'T-bin i the cross section {da/dET)i using 
the three sets of MRST99 with their correspond- 
ing as (A'/2)-values. By combining the as val- 
ues for the different E^ bins a value of = 
0.1224 ± 0.0002(stat) + 0.0022/ - 0.0019(exp) + 
0.0054/ - 0.0019(theor) was obtained. The major 
error sources are the jet energy scale (±1.5 % in 
as{Mz)) and the renormalization scale (-1-4.2/ — 
3.3% in as{Mz)); the uncertainties from the PDFs 
and hadronization are estimated to be less than 1%. 

Similarly, the dependence of on the energy 
scale was investigating by determining as from 
da/dET at different Et values. The resulting 
as{ET) values agree well with the predicted run- 
ning of as . 

Using the above data, ZEUS has measured the 



Figure 16: Mean subjet multiplicity as a function of 
ycut in inclusive jet production in DIS for Etj > 15 
GeV, > 125 GeV^ and -1< -qj < 2. 



dependence of the scaled invariant jet cross sections 
E^-Ejd^a/dpxjdpyjdpzj on the dimensionless vari- 
able XT — 2ETj /W.yp. The cross section measured 
at {W^p) of 180 and 255 GeV and averaged over 
-2 < ryJP < 0, agrees well with NLO QCD. In the 
quark parton model, the scaled cross section should 
be independent of W^p, whereas from QCD one ex- 
pects a violation of scaling due to the evolution of 
the structure functions and the running of as. As 
shown in fig. El the ratio of the scaled cross section 
at the two different energies W^p violates scaling as 
predicted by NLO QCD. 

3.2 Inclusive dijets 

Dijet production has been measured by ZEUS ^Sl 
and HI |16| under similar kinematical conditions as 
in single jet production. In both experiments the 
jets were reconstructed with the ILICA, requiring 
asymmetric cuts on Et in the lab frame. Both 
experiments cover a range of — 1 < 77^1^2 < 2.5 and 
0.1 < y < 0.9. The quantity x°'"' is introduced 

<^ = {ETjie-'^' + ETj2e~''')/2yEe (5) 

where y = W"^ /s, which in LO corresponds to the 
fraction of the photon momentum that contributes 
to the production of the two highest-i?T jets; cuts 
on allow the data sample to be enriched with 
direct processes (high x"^") or resolved processes 
(low x"*"*). The measurements are compared to 
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NLO QCD calculations with jin = jip equal set 
to the average i?T of the two outgoing partons. 
The PDF-parametrizations used are GRV-HO and 
AFG-HO for the photon and CTEQ5M1 for the 
proton; as{Mz) is set to 0.118. 
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Figure 18: Scaling violation in photoproduction (Q^ < 
1 GeV^): Ratio of measured scaled invariant jet cross 
sections (see text) for the two (W^p) shown and aver- 
aged over -2 < r^^ < 0. NLO QCD as in fig.lTTI 



in good agreement with the world average ffigs. I2ll 
and EH). 



Figure 17: Inclusive jet production in photoproduction 
(Q^ < 1 GeV^) integrated over — 1 < r/j < 2.5. LO and 
NLO: Klasen, Kleinwort, Kramer with — fip = Et, 
proton-PDFs: MRST99, photon-PDFs: GRV. 

The measured angular distributions (fig. El con- 
firm the steeper rise expected for resolved photon 
processes {x°'"' < 0.75) in comparison to direct 
photon processes {x°'"' > 0.75), indicating that the 
dynamics of the hard scattering process is reason- 
ably well described. For the HI data, the cross 
section as a function of (here identical to x"'"' ) , 
shown in fig. 1201 for two regions of E^^rnax, which 
represent two different factorization scales for the 
PDFs of the proton and photon, agrees with NLO 
predictions and varies only slightly with the dif- 
ferent parametrizations of the photon PDFs while 
ZEUS observes differences; it has to be noted, how- 
ever, that the different ii^T-cuts have been used (HI: 
Eti > 25 GeV, Et2 > 15 GeV; ZEUS: Eti > 14 
GeV, Et2>11 GeV). 

4 Summary of results on ag 

The results obtained from the recent measurements 
of jet and subjet cross sections in DIS and pho- 
toproduction on as{Mz) as well as on the scaling 
behaviour of as have achieved an accuracy compat- 
ible with results from other colliders. They are also 



5 Summary 

For large scales, and/or jet-i?T, the recent re- 
sults from HI and ZEUS on jet production in 
deep inelastic scattering and photoproduction are 
in good agreement with NLO QCD. In general, the 
hadronization corrections applied to the NLO re- 
sults are small and improve the description of the 
data only slightly. For decreasing scales and in spe- 
cific areas of phase space the corrections become 
large and are important in improving the descrip- 
tion of the data. 

In these studies, the longitudinally invariant 
kx cluster algorithm in its inclusive mode has 
proved reliable and, therefore, become the preferred 
method of jet identification in this field. The good 
agreement of as determinations performed with jet 
(and subjet) production at high and E^ with 
the world average supports this conclusion. 

The situation is worse in the small x / forward 
region. For x < 10^'* and small a breakdown 
of DGLAP is to be expected and indeed increas- 
ing discrepancies with DGLAP based predictions 
are observed even in NLO. In some cases in spite 
of small experimental errors the sizable theoretical 
uncertainties, however, do not yet allow safe con- 
clusions on a more satisfactory and comprehensive 
description of the data by modified pQCD models 
like BFKL or CCFM. 

Not only an improvement of experimental statis- 
tical and systematic errors, to come from HERA 
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Figure 19: Inclusive dijet production in photoproduc- 
tion {Q^ < 1 GeV^) as a function of the dijet angle 
0* in the parton-parton-cms for different cuts on the 
fractional momentum x°^^ of the photon participat- 
ing in the production of the two jets with the high- 
est transverse energy. Cuts: -ETji,2 > 14, 11 GeV; 
-1 < rij < 2.4; 134 < W-yp < 277 GeV. NLO: Frix- 
ione, Ridolfi with CTEQ5M for proton-PDFs and GRV 
resp. AFG for photon-PDFs. 



II, but also more accurate model calculations are 
needed to get closer to a satisfactory understand- 
ing of the interplay of soft and hard QCD processes 
in DIS and photoproduction. 
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Figure 22: Determination of as{Mz) and its scale de- 
pendence in recent HERA experiments. 
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